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ABSTRACT 

WIYN /HYDRA spectra in the Li 6708 A region have been obtained for 332 probable 
members of the old open cluster, NGC 6819. Preliminary analysis shows a pattern of 
Li depletion from the top of the turnoff to the base of the giant branch. Starting one 
magnitude below the level of the clump, all brighter giants have A(Li) below 1.0, with 
most having upper limits below 0.5. Star W007017, located below the first-ascent red 
giant bump is Li-rich with A(Li) = 2.3. As a highly probable single-star astrometric and 
radial-velocity cluster member, its discrepant asteroseismic membership could be a by- 
product of the processes that triggered Li-enhancement. Its color-magnitude diagram 
location is consistent with only one proposed enhanced mixing process among first- 
ascent red giants. 



Subject headings: open clusters and associations : individual (NGC 6819), stars : abun- 
dances, stars : chemically peculiar 
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Introduction 



Elemental abundances in stars supply insight into astrophysical processes governing Galactic 
chemical evolution and the structural evolution of stars. Few elements have garnered as much 
attention for their impact on cosmology and stellar evolution as Lithium. Created during Big Bang 
nucleosynthesis, the primordial Li abundance acts as a sensitive constraint for models of the early 
Universe. Easily destroyed in stellar interiors at T > 2.5 x 10 6 K, its variation with evolutionary 
phase is a diagnostic of the temporal significance of mixing and/or diffusion in stellar atmospheres. 
Not surprisingly, both areas suffer from di screpancies between standard theory and observations. 

for A(Li) (= 1 2 + log(N n /N h)) is three times higher 



The WMAP estimate (Larson et al. 2011 



than the observed limit in halo dwarfs ([Spite et al. 1120121 ) while disk field and open cluster stars 



show Li depletion where none should exist, e.g. the Li-dip among F dwarfs (jBoesgaard &; Tripicco 



1986) , or depletion ra tes that far exceed predicted values, e.g. among cooler dwarfs like the Sun 
(|xiong & Denglbood ). 



Assuming the WMAP-inferred value for the Big Bang Li abundance, a 4- fold increase in Li 
between Big Bang nucleosynthesis and solar system formation requires major Li sources on Galac- 
tic scales. Cosmic ray spallation remains the only viable source for production of 6 Li, but light 
element isoto pic abundances require a stellar origin for at least 70% of the Li associated with the 
solar system (jPrantzos 1 120121 . and references therein). The need to create Li in a non-destructive 
but observationally accessible environment has led to an emphasis on abundance anomalies among 
evolved stars, i.e. a search for the atmospheric signatures of the Cameron-Fowler (CF) mech- 
anism (ICameron &: Fowler Ill97ll ) generating Li in low-mass red giants via cool-bottom burning 
(jSackmann Boothroyd Ill999l) and in intermediate and higher mass asymptotic giant branch stars 
via hot-bottom burning (IKarakas 11201 ll . and references therein). To date, no consensus has emerged 



on the efficacy of either process. Searches for super-Li-rich giants, stars with A(Li) above the current 



Uttenthaler et al. 1120120 . typically have an incidence rate below 



up phase ( 


Kumar et al. 


2011 


( Smith 


20ld). In recent years 



giant s with anomalous Li abunda nces have been iden tified in the Gala ctic Bulge (ILebzelter et al. 



the field star popula tion (|Kumar et al. 



2012J) , dwarf spheroidal g alaxies (IKirbv et al. 1 12012| ) , t he thick disk ([Monaco et al. 



2011 



2011), and 



(|2OO0h and 



). While ICharbonnel Balachandran 
Kumar et al. I (|201ll ) have claimed that Li-rich giants concentrate in a narrow luminosity range 
associated with the red giant bump and/or the red giant clump, other investigations imply that 
Li-rich giants can be found almost anywhere along the red giant branch. This ambiguity makes 
assigning responsibility for Li enhancement difficult and is partially due to the range in mass, 
chemical composition, and evolutionary phase found among the mixed population of stars included 
in most surveys. Ideally, one would prefer only stars with a range in mass at a given age and 
composition or stars of a given mass and composition over a range in age, making star clusters 
ideal contexts for these studies. 



With the goal of probing stellar structure and evolution among low mass stars via atmo- 
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spheric Li, the authors have undertaken an extensive program to survey spectroscopically mem- 
bers of a key set of open clusters from the tip of the giant branch to as far down the main se- 
quence as technology allows. First result s have been published for the clusters NGC 3680 (age = 



1.75 Gvr) (lAnthonv-Twarog et al. 1120091 ) and NGC 6253 (3.0 Gyr) (|Anthonv-Twarog et al. Il20ld : 



Cummings et al. II2012I ). Among the clusters currently under analysis are NGC 7789 (1.5 Gyr) 



and NGC 6819 (2.3 Gyr), with over 300 stars observed in each. NGC 6819 has the distinction of 
lying wi thin the Kepler fiel d and remains the focus of asteroseismic studies reaching down the giant 



branch (jStello et al. 11201 ll ). The purpose of this Letter is to report the discovery of a unique Li-rich 
giant in NGC 6819. A detailed investigation of the theoretical implications of the Li pattern in 
NGC 6819, as well as the Li-rich giant, are bey ond the scope of this Letter and will be discussed in 
a future paper (jAnthonv-Twarog et al. II2013I ). We will, however, provide some context to clarify 
the unique nature of this highly probable cluster member. 



2. Observations and Reduction 



Spectra of 332 probable members of the open cluster NGC 6819 were obtained between 2010 
and 2012 using the WIYN 3.5m telecope at Kitt Peak National Observatory equipped with the 
HYDRA multi-object fiber-based spectrograph. Six configurations were used to extend the sample 
from V = 11 to 16.5. Complete details of the observation and reduction procedures will be pro- 
vided as part of a comp r ehens ive discussion of the spect ra and Li abundan c es for all 332 stars in 



Anthony- Twarog et al. I (|2013l ) . The reader is referred to ICummings et al. I ([2012]) for insight into 



the calibration and reduction steps used within the project. 

The cluster spectra cover a wavelength range from 6550 to 6840 A with a resolution of 0.5 A 
= 2.7 pixels. Due to the plethora of lines for cool giants at the metallicity of NGC 6819, we have 
estimated the average S/N per pixel in this wavelength region based on the accumulated photon 
statistics. While these probably optimistic S/N ratios per pixel for the entire sample range from 
below 70 to several hundred, the evolved stars discussed in the current Letter have S/N per pixel 
above 100, with a majority above 200. 



3. The Color-Magnitude Diagram and Li Evolution 



Sam ple selection was ba sed on the comprehensive photometric and radial- velocity study of the 
cluster by lHole et al. I (J2009j)(H09). All stars chosen were classed by H09 as radial- velocity members 
with a probability above 50%. The color-magnitude diagram (CMD) using photometry from Table 
2 of H09 for our 332 stars is shown in Fig. 1. Th e magenta starred point is the proposed Li-rich 
giant. Superposed is a Y 2 ( Demarque et al. 1120041 ) isochrone with an age of 2.25 Gyr and [Fe/H] = 
+0.09, adjusted for E(B — V) = 0.14 and (to — M ) = 12.5. The Y 2 isochrones have been adopted 



for consistency with previous cluster studies (jAnthonv-Twarog et al. 1120091 . 



2010; 



Carraro et al. 
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2011). 



We have adopted the reddening value and metallicity derived by iBragaglia et al. I (|200ll ) from 
high-dispersion spectroscopy of four clump giants. E(B — V) = 0.14 falls within the range E(B — V) 
of 0.12 to 0.16 found in published studi es of the cluster using method s such as differential cluster 
CMD match es to isochrones and M 67 (jRosvick fe VandenBerg 1 1 19981 } . Washington p hotoelectric 



phot o metry (ICanterna et al. I [198(f) , and UBV photoelectric and photographic data (jBurkhead 
197lh. iKalirai et al. I ( 200ll) adopt E(B - V) = 0.10. The adopted metallicity ([Fe/H]= 0.09) from 



Bragaglia et al. I (120011) is cons i stent w ith the adopted reddening; moderate-dispersion spectroscopy 



of seven giants by iFriel et al. I (|2002i ) usin g E(B — V) = 0.16 gives [Fe/H] = 0.04 if their scale is 



adjusted so that M67 has [Fe/H] = 0.0 (jTwarog et al. 1119971 ). Ultimately, our conclusions are 



remarkably insensitive to the value of the reddening (see Sec. 4). 

With the reddening and metallicity known, T e f f was derived for e ach star using the B— V of Fig 



1 and the color-temperature- metallicity relations for cool giants from lAlonso et al. I (|1999l ). Surface 



gravities were obtained by mapping V to log g values established by the isochrone com parison (Fig 
1). For 21 stars with mass and radius information from the Kepler analyses by IStello et al. 



(|201ll ). we compared the isochrone-derived surface gravities with those from asteroseismology and 
accordingly incremented the isochrone-implied values by 0.08. 

Equivalent (EW) measures were made by multiple analysts using SPLOT to measure lines 
within 10 A of the Li 6708 A feature, then combined to a common measurement system. After 
subtracting the expected contr ibution to the 6708 A feature that is due to Fe I at 6707.45 A 
(jSteinhauer Delivannis II2004I ). we estimate the Li ab undance by interpolating within a grid of 
temperature and Li abundances from ISteinhauer I (|2003l ). 



While this approach works well for stars at the turnoff, it is less satisfactory for evolved 
stars whose temperatures are near the low temperature limit of the grid's applicability. Also, the 
line-dense spectra of the giants complicates placement of the stellar continuum near the typically 
weak Li line, so we utilized spectrum synthesis techniques. An appropriate model atmosphere for 
each star is specified by its color-based T e ff, its U-dependent log g an d a microturbulent vel ocity 



estimated from surface gravity as follows: vt = 2.0 — 0.2 log g following ICarreta et al. I (|2004l ) and 



McWilliam fe Bernstein I (|2008l ). Using MOOG, comparison of the synthesized spectrum, for which 
the Li abundance may be varied interactively, to the observed spectrum leads to an optimal A(Li). 
We stress that for the majority of giants, the 6708 A feature is dominated by the nearby Fe I line, 
so only an upper limit is possible for Li. The upper limits on A(Li) are conservative in that the 
limit would need to be reduced by at least 0.3 dex for the synthetic spectrum to approach the 
observed line. 

Fig. 2a shows the CMD for all evolved stars redder than B — V = 0.70. The caption summarizes 
symbol keys that provide information on binarity as well as highlighting asteroseismic non- members, 
subgiants and the proposed Li-rich giant. Fig. 2b shows the same stars with the measured A(Li) 
as the primary axis. The mean value for the stars at the turnoff above the Li-dip is marked by the 
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vertical arrow at A(Li) = 3.3. Fig. 3 illustrates the spectra for W007017, with an upper limit of 
Vsini = 8.3 ± 0.3 km/sec from IRAF's fxcor, and a cluster member with nearly identical color and 
temperature, W005011. 



4. Discussion and Conclusions 



Ignoring binaries and/or potential non-members, the evolutionary patte rn in Fig. 2 resembles 



that outlined in the younger clusters NGC 752, NGC 3680, and IC 4651 (|Pasauini et al. 112004 



Anthonv-Twarog et al. 1 12009 ) . Unevolved stars at the turnoff on the hot side of the Li-dip region 
have the solar system value of A(Li) = 3.3 and are thus inferred to be depleted minimally, if at 
all. As stars evolve across the subgiant bra nch, the SCZ deepens drastically. A(Li) is depleted 
by standard dilution (jDelivannis et al. Ill990l ) and/or non-standard mixing to below 1.5 when the 
SCZ reaches its maximum depth. For stars with V < 14, there is a precipitous drop in A(Li); 
typical upper limits of A(Li) = 0.5 or less are found. As expected due to the T e ff sensitivity 
of the Li feature (and higher S/N), the upper limits implied by a given line strength decline as 
stars become redder. Only 6 stars with V < 14, including two known binaries, have measurable 
A(Li). All but one have A(Li) below 1.0. The Li-rich giant, W007017 (WOCS ID from H09), has 
A(Li) = 2.3 and is located below the red giant clump and the red giant bump by ~0.5 mag. From 
standard stellar evolution, this implies that the star must be a first-ascent red giant with an inert 
He-core surrounded by a H-burning shell. The latter comparison is critical because the red giant 
bump marks the phase where the H-burning shell crosses the chemical discontinuity left behind 
by the bottom of the convective envelope at the end of the first dredge-up. Above this luminosity 
level, lower mass stars are believed to undergo extra mixing in their convecti yely stable radiative 



zones between the H-burning shell and the bottom of the convective envelope (jRvan &: Delivannis 



19951 ). leading to enhanced atmospheric Li as a by-product of the CF mechanism. The CMD 
position of this star is incompa tible with all current models for Li enhancement in low mass giants 
(jCharbonnel &; Lagarde II201QJ . and references therein) except one. Before discussing the exception, 
a few issues relevant to the paradoxical CMD location of the Li-rich giant should be addressed. 



a) Is W007017 a member of N GC 6819? It is a p roper motion (90%) ([Sanders Ill972l ) member, 
a single-star radial-velocity (95%) ( Hole et al. Il2009l )member. and a photometric member (located 



at the red edge of the RGB). The first two together imply near-certainty of cluster membership; 
furthermore, it is even m ore improbab l e that a non-member would meet all three of these criteria 
by coincidence. And yet, IStello et al. I (120111 ) have classified it as a seismic non-member. Cluster 
membership based upon seismic oscillations is a measure of the degree of similarity in the oscillation 
pattern among stars at the same position/evolutionary phase in the CMD. Four of the six seismic 
non-members in NGC 6819 are also photometric non-members. The two photometric members, 
W007021 and W007017, lie along the red edge of the first-ascent giant branch, 0.5 to 0.7 mag below 
the level of the clump. Of the six seismic non-members, W007021 exhibits the greatest deviation 
from the relations defined by the cluster members while W007017 shows the smallest. Seismic 
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oscillations have the potential to identify as non-members any field interlopers whose kinematics 
mimic that of the c luster, but also any exotic stars, such as remnants of strong dynamic interactions 



(jStello et al. 11201 ll ) because standard evolution is assumed for cluster stars when estimating seismic 
parameters. While it isn't possible to entirely rule out the possibility that W007017 is be a Li-rich 
field star mimicking the cluster space motion and photometric sequence, it seems far more probable 
that its seismic deviation from the cluster relations reflects the structural mechanism that led to 
its anomalously high atmospheric Li abundance. 

b) How reliable are the cluster /stellar parameters? A(Li) depends upon the derived T e ff, log 
g, and vt, in turn dependent upon the adopted [Fe/H], reddening, distance modulus, and age. To 
test the sensitivity of our conclusions to these parameters, A(Li) was rederived using synthetic 
spectra with T e ff altered by ±100 K, log g by ±0.1, and [Fe/H] by ±0.1 dex. While systematic 
shifts occurred in the absolute scale of Fig. 2b, the relative A (Li) ranking of the stars remained the 
same within the measurement errors. As an example, lowering T e ff by 100 K on average dropped 
A(Li) by 0.1 dex, with a smaller shift of 0.05 dex near the turnoff, increasing to 0.15 dex among 
the cooler giants; W007017 still remained anomalously rich (A(Li) = 2.2) compared to the evolved 
stars brighter than V = 14.0. A related question is tied to the discrepancy between the position of 
W007017 in the CMD and the red giant bump. As the cluster ages and lower mass stars evolve up 
the giant bra nch, the luminosity of th e bump location drops and the luminosity range of the bump 



grows larger (jDemarque et al. 1 12004). Would increasing the age of the cluster place the Li-rich 
giant above the bump? As implied by Fig. 2a, increasing the age of the cluster from 2.25 Gyr to 
2.5 Gyr would require a reddening estimate reduced to E(B — V) = 0.10 to match the turnoff color. 
In turn, lowering th e reddening would lo wer the apparent modulus to 12.3, the same combination 



of values derived by iKalirai et al. I (120011 ) using fits to the Hyades cluster. So, while the My limit 
for the bump would be about 0.2 mag fainter for the older giant branch, the V limit would remain 
near 13.0, unchanged from the current estimate. 

c) Is W007017 a composite system? Of the 6 stars with measurable Li, 5 lie within 1 mag 
of each other near the clump and, of these, two are known binaries. While the measurements to 
date number 5, the radial- velocity data of H09 and the current study show no evidence for radial- 
velocity variations or unusual line broadening, so the data are consistent with a single star. It 
should be noted that two fainter, bluer stars lie within 5" of W007017. The brighter star at V = 
16.69 and B - V = 0.81 is 3.3" away while the fainter star at V = 17.18 and B — V = 0.70 is 4.9" 
removed, while W007017 has V = 13.57 and B — V = 1.24. Since the CCD photometry used in the 
determination of V and B — V for W007017 is based upon point-spread-function fitting techniques, 
it is unlikely that the colors for W007017 are seriously impacted by the presence of this pair. If, 
however, all three stars were erroneously measured as one system, the composite would be brighter 
by 0.1 mag and bluer by 0.05. More realistic mixtures of the contaminating stars scaled by the 
square of their distance from W007017 lead to upper limits on the impact of these two stars of 
-0.05 mag in V and -0.03 mag in B — V. Therefore, in the limit, W007017 could have true colors 
closer to V = 13.62 and B — V = 1.27. A shift in B - V by 0.03 alters the T eff of W007017 by 66 
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K and lowers the derived A(Li) by less than 0.1 dex. For the seismic analysis, the c ont aminat ing 
pair would increase the noise of the signal but not alter the oscillation frequencies (jStello et al 



201 ih . 



Although the evidence for cluster membership is compelling, W007017 is in any case a member 
of the rare class of red giants with anomalously high atmospheric Li. If a member of NGC 6819, its 
position in the CMD places it below both the red giant clump and red giant bump for first-ascent 
red giants, leaving the obvious question of the origin of the enhanced Li. If t here is no interactin g 
companion, an extremely salient possibility emerges from the recent work of iDenissenkov I (|2012l ). 
In attempting to exp l ain the appar e nt di stribution of red giants below the bump as compiled by 



Kumar et al. I (|201ll ) , IDenissenkov I ([2012]) found that single star models undergoing fast internal 
rotation could experience extra turbulent mixing, resulting in enhanced surface Li at the level found 
in W007017. There could also be extra heat transport, affecting the structure in the radiative zone 
above the H-burning shell (by increasing dlnT/dlnP) as the star evolves brighter than the red giant 
bump luminosity, causing the star to zigzag toward lower luminosity along a cooler evolutionary 
track before reascending the giant branch. Thus, Li-rich giants on the first-ascent of the giant branch 
should be positioned redward of the standard giant branch, fainter than the red giant bump, and 
should have an atmospheric structure that differs significantly from standard giants ascending to 
the red giant bump for the first time. Such stars could be relatively rare, since only the most 
rapidly rotating stars would become Li-rich and execute zigzags. Whether this scenario applies to 
W007017 remains to be seen, but the unique properties of this star indicate that additional scrutiny, 
particularly searching for spectroscopic signatures of exceptional mixing, such as, for example, the 
12 C/ 13 C ratio, could be beneficial to our understanding of this rare class of objects. 
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support provided by then grad students Ryan Maderak and Jeff Cummings of IU in collecting the 
spectra. NSF support for this project was provided to ER through NSF grant AST-0850564 via 
the CSUURE REU program supervised by Eric Sandquist, to BJAT and BAT through NSF grant 
AST-1211621, and to CPD through NSF grant AST-1211699. Useful feedback on the interpreta- 
tion of the asteroseismology was supplied by Ron Gilliland and Dennis Stello. Extensive use was 
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Fig. 1. — CMD of stars observed spectroscopically wi thin NGC 6819. Open circles are single- 
star members, as verified by t he radial-velocity s tudy by lHole et al. I (J2009J). Filled symbols denote 
probable binaries, as noted by Hole et al. (j2009l ). black and the present study, red . The starred ma- 
genta point marks the position of the proposed Li-rich giant. Superposed is a Y 2 (jDemarque et al. 
20041 ) isochrone with an age of 2.25 Gyr and [Fe/H] = +0.09, adjusted for E(B - V) = 0.14 and 
(m-M) = 12.5. 
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Fig. 2.— Left panel: CMD of the evolved stars ((S - V) > 0.70) with derived A(Li). Open 
circles are single-star members with measurable Li; open triangles are stars with upper limits for 
A(Li). Filled symbols are probable binaries fro H09 (bla c k) an d the current study (red). Green 
symbols are asteroseismic non- members from lStello et al. I (|201ll ). Blue symbols delineate subgiant 
stars. The magenta starred point is the Li-rich giant. Right panel: A(Li) for stars along the CMD. 
Symbols have the same meaning as in the left panel. 
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Fig. 3.— Spectra near the 6707.8 A line for red giants W007017 (blue) and W005011. Note the Fe 
I line at 6707.45 A in W005011. 



